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^ : Abstract 

O 
CD 

The approach of multi-dimensional SUSY Quantum Mechanics is used in an explicit con- 
struction of exactly solvable 3-body ( and quasi-exactly-solvable iV-body ) matrix problems on 
| a line. From intertwining relations with time-dependent operators, we build exactly solvable 

non-stationary scalar and 2x2 matrix 3-body models which are time-dependent extensions of 
the Calogero model. Finally, we investigate the invariant operators associated to these systems. 

t> 
(N 

O ■ !• Introduction 
O 

During the last three decades exactly solvable iV-body problems have provided useful 
tools to investigate formal algebraic properties with applications to different branches of 
Physics. The most widely studied model is the so called Calogero model |l| and its var- 
Q-f ious generalizations which essentially are many-body extensions of the one-dimensional 
singular harmomic oscillator model. Calogero-like models have been developed incor- 



porating different root systems 0, g-deformations [Q], PT-symmetric generalizations ||, 
many-body forces ||, multi-dimensions || and internal degrees of freedom with potentials 
which couple them (matrix potentials) J7J. Even if coupled channel problems in general 
have a continuum || and discrete spectrum, in the Calogero-like models (with harmonic 
attraction) the spectrum is purely discrete reflecting an essentially confining dynamics. 
Physical applications of this dynamics have been elaborated in the context of localized 
systems like, in particular, Paul traps and quantum dots Q. 

An important generalization of the Calogero model is its supersymmetrization JIQ| . 
Supersymmetric Quantum Mechanics (SUSY QM) Jll| is a suitable framework to dis- 



cover and investigate one- dimensional, one-particle exactly-solvable models [T2|]. The 
same strategy applies to systems with multiple degrees of freedom. Multi-dimensional 
SUSY QM, first constructed in fl3|], leads to a superHamiltonian which includes a chain 



of matrix Hamiltonians (cf. coupled channels or internal degrees of freedom, like spin 
fl4|). Supersymmetry ensures that the spectral properties and the eigenfunctions of the 
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Hamiltonians belonging to the chain are algebraically interrelated. This analysis can be 
reinterpreted with reference to a one-dimensional multi-particle problem enlarging classes 
of many-body exactly solvable problems in a similar way as for the one- dimensional one- 
particle problems [pl| . 

We start from the superpotential of the Calogero system which corresponds to two 
exactly-solvable scalar components of the superHamiltonian. They are intertwined to the 
neighbouring matrix components of the superHamiltonian. This implies that a part of 
the spectrum for both matrix potentials and the corresponding wave functions are known. 
Thus from solvable scalar models by supersymmetric techniques quasi-exactly-solvable 
matrix problems are generated. This approach generates matrix iV-particle models which 
can be considered in the context of recently constructed scalar quasi-exactly-solvable [16 
and so called partially solvable [IT]] models. 



In Section 2 we review the basic aspects of multi-dimensional SUSY QM |13[ and 
introduce its reinterpretation in terms of multi-particle one-dimensional SUSY QM (see 
details in the recent paper |15[]). In particular, we focus attention on models with exactly- 
solvable scalar components of the superHamiltonian. 

In Section 3 starting from the Calogero model |l| we analyze its SUSY extension 
which includes quasi-exactly-solvable |TB[ iV— particle matrix models. Furthermore, we 



study 3-body problems in detail because the properties of the chain of the components 
of the superHamiltonian simplify considerably so that the spectrum and the wave func- 
tions for the (only) matrix Hamiltonian are fully determined from those of (two) scalar 
Hamiltonians. 

In addition to exactly- solvable stationary problems we consider also time-dependent 
potentials and, correspondingly, exactly-solvable time-dependent problems. In the context 
of one- dimensional one-particle SUSY QM (and Darboux transformation) such problems 
were investigated in [19| . In non-stationary SUSY QM supercharges (of first and second 
order in space derivatives) commute with the non-stationary Schrodinger superoperator 
and intertwine consecutive components of the supersymmetric chain. Following meth- 
ods developed in recent investigations of the time-dependent harmonic oscillator model 
and its generalizations f2"0|| , in Section 4 we construct time-dependent 3-particle solvable 
problems. In this Section we achieve our main goal after having prepared in the previ- 
ous Sections the relevant framework. These results can be interpreted as time-dependent 
generalization of the SUSY Calogero model, which can be shown to be solvable by intro- 
ducing unitary intertwining operators (non-polynomial in derivatives). An extension of 
this method to the iV-body Calogero model described in Section 2 leads to time-dependent 
quasi-exactly-solvable matrix models. While the Calogero model is a many-body gener- 
alization of the "singular" harmonic oscillator model, its time-dependent version, which 
we study, are correspondingly extensions of the oscillator problem with time- dependent 
parameters. This last problem has attracted much interest in the literature pD |, [21] and 
has applications in different areas of Physics. 
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2. Multi-dimensional SUSY QM and iV-particle quasi- exactly- 
solvable stationary problems 

The supersymmetric quantum system for arbitrary number of dimensions N consists 
]T3| of the superHamiltonian H s and the supercharges Q ± with the algebra (SUSY QM 
algebra) : 

H S = {Q + ,Q~}, (1) 
(Q+f = (g-)^ = o, (2) 

[H s ,Q ± ] = 0. (3) 

An explicit realization is given byfj 

i N N N 

H s = - (-A + ]T(cW) 2 - AW ) + E <■■; CjihiJjW: A = J2 didi, d t = d/dx i; (4) 
— i=i *)i=i i=i 



1 N 

Q ± = ^E#(±^ + W ) ( 5 ) 



where ^j, ■?/>/" are standard fermionic operators: 

{Vi, ^} = 0, {^ + , V/} = 0, = <5y, (6) 

The dynamics of a particular SUSY QM system is determined by a superpotential W, 
which depends on iV coordinates (x\, . . . , x^). 

In general, solvable scalar models in multi-dimensional Quantum Mechanics for one 
particle admit simple separation of variables and are, therefore, reducible to one-dimensional 
problems. For this reason from now on we will alternatively interpret multi-dimensional 
SUSY QM as a multi-particle problem on a line because one knows classes of solvable 
models (Calogero model, Sutherland model and others [[|) which do not admit such a 
straightforward separation. 

For iV-particle systems on a line it is natural to consider [15] superpotentials with a 
separable centre-of-mass motion (CMM), satisfying the condition: 

( \ N 

W(xi, ...,x N ) = w(xi, ...,x N ) + W {x 1 + ...+x N j; 22 d J w ( x ^ ■■■,x N ) =0, (7) 

i.e. the first term w(xi, . . . , x^) does not depend on J^iLi x i- We will restrict ourselves to 
superpotentials (0) with Wq = 0, or, equivalently, J2k=i dkW(x±, ...,xn) = 0. 

For the superpotentials (|7p one can use the well-known Jacobi coordinates^ |22| : 

Vb = , (xi + ...+x b - bx b+1 ); (8) 

'b(b + 1) 

N 

X 



i=l 



^^Here and below the indices k, . . . range from 1 to N. 

Prom this moment on, the variables denoted by letters a, b, c, . . . range from 1 to (N — 1). 
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or, shortly, yu = J2i=iRki x h where the orthogonal matrix R is determined by @. For 
the supersymmetric systems we introduce also the fermionic analogues^ of the Jacobi 
variables: 

(V>i + . . . + ^6 - bipb+i); 



6(6 + 1 



1 N 



which satisfy the canonical anticommutation relations: 

{0 fc , 0,1 = 0, {0+0+} = 0, {0 fc ,0+}=^. (9) 
In terms of the Jacobi variables the supercharges (|5|) can be rewritten as: 

Because 

{5 ± ,^/-} = 0, (10) 
the free motion of center-of-mass in the superHamiltonian can be separated: 

id 2 

where 

l N - 1 / d 2 /dw\ 2 d 2 w\ N_1 d 2 w 
h = {<?+, q } = - E - 0-2 + «~ " ITT + E tifc-z-z- 



fl5 = {0 + ,0 _ } = /*-sCT. ( n ) 



(12) 



is (iV — l)-dimensional superHamiltonian expressed in Jacobi variables yi, yn-i- in the 
following we will consider only this reduced superHamiltonian h. 
The operator h acting in the fermionic Fock space 



|0 >; M < N; 6< < for i < j, (13) 



generated by fermionic creation operators 0+, conserves the corresponding fermionic num- 
ber 

N-l 

[h, N F ] = with N F =J2 ttK ( 14 ) 

6=1 



3 The use of these variables has been instrumental |L5j to clarify the role of the Permutation Group Sn 
in SUSY QM. 
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Therefore, in the basis (|T3|) it has |13j a block-diagonal form: 



h 



















h& . . . 
















h (N-2) 







V o 


0.. 







V 



(15) 



where the matrix operators of dimension Cf}_ 1 x C^}_ 1 are the projections of h onto 
the subspaces with fixed fermionic number Np = M. These components are standard 
Schrodinger operators with matrix potentials and can be obtained from fll2|) by a suitable 
matrix realization of the fermionic variables 

The supercharge q + increases the fermionic number from M to M + 1 and has the 
under-diagonal structure: 





( ° 







9(0,1) 





q + = 





9(1,2) 




V o 





















(16) 



Q(N-2,N-1) 0/ 



Similarly, q = (q + y is an over-diagonal matrix operator with nonzero elements 

Q(M+i,M) = ^9(Af,M+l) 

Superinvariance (|j) of the superHamiltonian corresponds, in components, to the in- 
tertwining relations: 

hq + = q+h * ^ (M+1) ? ( V + i) = *lM,M + i) hiM) 
q-h = hq- & q^ M+ x,M) h{M+1) = h{M)( l{M+i,M) 

These relations lead [|13| to important connections between spectra and eigenfunctions 
of "neighbouring" Hamiltonians, with fermionic numbers differing by 1. In particular, 
Q(m m+i) ma P s eigenfunctions of onto those of /i( M+1 ) with the same energy E K : 



Q(M,M+1) 



(17) 



Analogously, q^ M M _ 1 ) maps eigenfunctions of h^ M ^ onto those of with the same 
value of energy (see details in [|L3|j ). 



In particular, the spectrum of the matrix (N — 1) x (N — 1) Hamiltonian consists 
of two portions, one of which coincides with the spectrum of the scalar Hamiltonian . 
Thus if the scalar problem with is solvable the matrix problem with h^j) becomes 
quasi-exactly solvable Hl8| . Similarly, the matrix Hamiltonian hS N ~ 2 ^ is also quasi-exactly 
solvable provided the last (scalar) Hamiltonian h^ N ~ lS) is exactly-solvable. 
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3. Stationary solutions of 3-body problem with internal degrees 
of freedom 



As a realization of what we presented in Section 2, we provide an explicit construction 
for the iV-body Calogero model. Substituting the superpotential [] which depends only 
on first (N — 1) bosonic Jacobi coordinates yi,y 2 , ■■■,Vn-i' 

N 7 N 

W{x 1 ,x 2 ,...,x N ) = a J2( x i- X j) +2 E xj\ =w{y 1 ,y 2 ,...,y N -i)- (18) 

into @, after some manipulations we obtain apart from a constant energy shift: 

1 N 1 N 1 

Hs = -^+4a 2 N E (^-^) 2 + ^7(7 + 1) E , _ ~y 2 + 

N 1 iV 



jv 11 

+ 7 E ^ViT— la-T E ^T^r^- (is) 



After subtraction of the free center-of-mass motion (|TT| ) one obtains | a reduced Hamilto- 
nian h from (0), with the superpotential w(yi, y 2 , 2/jv-i)- The expression for scalar 
can be derived from the superHamiltonian (|l^) by taking into account that the fermionic 
terms vanish in the subspace with Np = : 

1 N 1 N 1 

^ (0) = -~ A i N ~ 1] + ^ N E 0* - + ^7(7 + 1) E T-T^i- (20) 

It corresponds to the well-known exactly solvable iV-body Calogero model 0,0. As was 
discussed in the end of the previous Section, the matrix Hamiltonian is thus quasi- 
exactly solvable and the associated part of its energy levels coincides with oscillator-like 
spectrum of (p0|). 

The last scalar component /i^ -1 ) of the superHamiltonian ( fT9|) is obtained by its 
reduction to the subspace of (|T^) with maximal fermionic occupation number Np = 
(N — 1). Only the last fermionic term in fll9|) is effective and /i^ -1 ) coincides^ with 
after the 7 into (—7) replacement [J. It is clear that exactly solvability of h^' 1 ^ leads 
again to quasi-exactly solvability of the matrix Hamiltonian h^ N ~ 2 '. 

For N = 4 the chain of (p~5|) consists of two scalar Calogero Hamiltonians hS Q \ 



and two matrix 3x3 Hamiltonians and hf^, where for example |L3 



hft = + ^{d % wf + i I ' 2dld 2 w {d 2 2 - d\ - di)w 20- 2( ),,r I (21 ) 



(df — d 2 — d%)w 2did 2 w 2<9i<9 3 u> 



2d!d 3 w 2d 2 d 3 w (Sf - d\ - d\ 



4 This form for W is suggested by the ground state wave function of the conventional Calogero model 
(see (p0[)) and it is known to be a particular choice among possible alternatives. 
5 From now on we will use the notation di for d/dyt. 

6 Let us note that the eigenfunctions of h^ ' and h^ -1 ' are not connected directly by supercharges q^, 
contrary to the hypothesis of the paper p3| in the context of Calogero-like models. In this connection, it 
was recently found |24[ that their eigenfunctions are related by a Dunkl-like differential operator. 

7 Note that the ground state energy of the Calogero model depends on 7. 
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and h ik has a similar structure. The Hamiltonian (|2ip is intertwined to /i 1 - - 1 by 1 - ) = 
(Ajf, Aj, A% ), where = (A+)t = 775 (9i + ^^(2/1, 2/2, 2/3))- Therefore since /i < -°' ) is solv- 



able, is quasi-exactly solvable. Similar considerations hold concerning the intertwining 



of and h®. The " non-quasi-exactly solvable" portions of and hff} coincide [|13[ 
because of an additional intertwining between them. 

It is clear that when the matrix operator happens to coalesce with hS N ~ 2 \ the 
quasi-exactly-solvable matrix problem becomes exactly solvable. This is the case for 
N = 3 Calogero model. We now consider the standard Calogero Hamiltonian for three 
particles on a line with repulsive singular terms. In terms of Jacobi coordinates 



Vi 



xi - x 2 



2/2 



X\ + x 2 - 2x 3 



V2 ' yz Vq 

the superpotential w(yi,y 2 ) up to an irrelevant constant has the form: 



w{yi,y 2 ) = Qa(y l + y 2 ) +7 In {y^-yi + — y 2 )(--yi + — y 2 )\ 



(22) 



The Hamiltonian 



can be rewritten as: 



h {0) = -^(9 1 2 + a 2 2 ) + i 7 ( 7 + l){l4 



vl 



+ 72a 2 (y 2 1 +y 2 2 ), 
or, equivalently: 



+ 



fi °) _ A^A X + A\A 2 , 
where (Ai , A 2 ) are the components of the vector operator 

9(o,i) = (A > A 2 ), 
which can be expressed in terms of the superpotential as: 



} + 



(23) 
(24) 

(25) 



Al 



-j=(d 1 + diwiy^y^) 



+ 12^ + 1 - 



2 

— + 



(A+)t = (d 2 + d 2 w{y u y 2 )) 



V2 
V2\ 



d 2 + I2ay 2 + 



+ 



^2 



The Hamiltonian is not symmetric under the exchange of variables y±,y 2 , how- 
ever its wave functions can be obtained from the well known wave functions of Calogero 
Hamiltonian [J], which are symmetric under the permutations of x^ (i — 1, 2, 3). 
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According to Section 2, the Hamiltonian generates a chain which includes a second 
scalar Hamiltonian defined (apart from a constant) by: 

/»« = B +B^ +B+B 2 = - l -{d\ + di) + 

+ 57(7-i){A + 7i — V^ + r^ — ^r^> + 72a2 ^i + ^) ; (26) 

where we have introduced the operators 5^ = eikA^; ei2 = — £21 = 1; en = £22 = 0. 
Also included in the chain is the 2x2 matrix Hamiltonian: 

,(*) - /i- /i + 



ftp) 



i (M+ i [(arf -^] + i(^ 



02 W dflU 



did 2 w 



l -{d\ + «9 2 2 ) + 72a 2 {y{ + y 2 2 ) + 36a 7 + ■ 



2 S ~ J 1 " ' " " ' " 1 </r 



where cij are Pauli matrices. 

The above Hamiltonians and are intertwined with h^ 1 ' : 

fi^Af = Afrhki; A l h^ > —h lk A k \ 

h®B+ = B+h$; Brh®=h£ ) Bt ] l,k = 1,2, (28) 

This chain of Hamiltonians h^°\ h\ k , determines the superHamiltonian as a Schrodinger- 
like operator with 4x4 matrix potential of block-diagonal form. Intertwining relations 
( p8|) lead to interrelations between spectra and eigenfunctions of the chain Hamiltonians. 
Apart from possible zero modes of Af , B^, the spectrum of 2 x 2 matrix Hamiltonian 
is formed by two parts, coinciding with the spectra of the scalar Hamiltonians and 
h^ 2 \ correspondingly. Their eigenfunctionsf] are connected by the intertwining operators: 

*£\e®) ~ A^¥°\E^y, v£\eW) ~ B k ¥ 2 \E^y, 

*<f>)(E<®) ~ A+*£\eW); ^ 2 \E^) ~ B+*£\EW). (29) 

Thus all (up to zero modes of the Af, Bf) eigenvectors of the matrix Hamiltonian hS 1 ' 
are expressed in terms of the Calogero wave functions. 

In summary, we have used the framework of SUSY QM in order to derive an exactly 
solvable 2x2 matrix model, the spectrum of which is divided into two parts, each one 
coinciding with the spectrum of a scalar Calogero Hamiltonian. The reason why the spec- 
trum of the matrix model (£27|) is still completely discrete can be found in the dominance 
of the confining scalar interaction over the coupling of internal degrees of freedom which 
is asymptotically decreasing. This matrix problem in a non-trivial way is related to a 
system of independent harmonic oscillators ||25|| , but is not diagonalizable by standard 
transformations like rotations. 

8^(1) (£(0)) 

are the components (I = 1,2) of two-component vector eigenfunctions of the matrix 
Hamiltonian M 1 ). 
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4. Time-dependent exactly solvable 3-body matrix problems 

In this Section we will achieve the goal of obtaining scalar and matrix time-dependent 
exactly(quasi-exactly) solvable models and invariant operators. We start from a gen- 
eral time-dependent intertwining relations which connect two time- dependent Schrodinger 
equations (TDSE) 0, one of them with a time-independent exactly-solvable Hamiltonian. 
If H(y) is an exactly solvable Hamiltonian and H(y) , ip n (y) = E n ip n (y), the intertwining 
relation with a known operator U(y, t) : 

(id t - H(y,t))U(y,t) = U(y,t)(idt ~ H{$)) (30) 

leads to an exactly solvable time-dependent problem. All the solutions of 

(id t -H(y,t))$(y,t)=0 

can be written as U(y,t)^(y,t), where ^f(y,t) = J2'^=o c n e ~ lEnt 4'n{y) is a generic time- 
dependent solution of equation {id t — H(y))^/(y,t) = 0. 



For the one-dimensional case intertwining relations ([3(]) were investigated in (19| for 
differential operators U(y, t) of first and second order in derivatives. While in the one- 
dimensional problem a wide class of solutions was found, a straightforward extension to 
the two-dimensional case does not appear to be obvious. In this case it is more effective to 
study operators U (y, t) which can be written as products of two unitary pseudo-differential 
(of infinite order in derivatives) operators of the form 



U(y,t)=exp{ia(t)J2y-}-exp{b(t)J2(y l d t + d t y t )}, (31) 

i i 

where a(t), b(t)- are arbitrary external time- dependent real functions. These operators 
have no zero modes. The intertwining relation (p0|) leads to: 



H(y,t)) = U(y,t)H(y)U-\y,t) +l (^^-)U- l (y,t). (32) 

In the supersymmetric framework (Sections 2 and 3) for each Hamiltonian of the chain 
one can choose the real valued coefficient functions a^ M '(t), b^ M \t) independently for the 
different values of M. Under these unitary transformations the Jacobi canonical 

variables transform as: 

Vi - f/W^(^ M ))- 1 =y J -ex P {26W(t)}; 
Pi = -idi - U^ Pl (U (M ^r 1 = (p t -2a^ M \t)y i )-exp{-2b^ M \t)}, 

and the so called gauge term in (|3"2"D reads: 



i i 

(33) 



9 Both Hamiltonians are assumed to be hermitian. 
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After setting up the general framework of time-dependent intertwining of TDSE, we 
apply it to the Hamiltonians hS M "> of the Calogero superchain of Section 3. In particular, 
we identify H(y) with the elements of the 3-body Calogero chain M = 0, 1, 2 and 
generate a time-dependent chain. In general, the time- dependent Hamiltonians acquire 
new terms linear in momenta and time dependent coefficients in all terms: 



-e 
2 



-46(°) (t) 



i=l,2 



E Pi ~ (« (0) (t)e- 4b(0) « + W(t) ) £ (y iPi +p iyi ) + 



i=l,2 



+ 



2(a(°)(t))V 4b(0) « + 72aV 6(0) « + Aa^(t)b^(t) - a^(t)) E vl + 



i=l,2 



+ ie-^ 0) W 7 (7 + l) 



y 2 



(34) 



The second scalar Hamiltonian (y, t) results from (|2~6| ) with a similar construction. 
The matrix Hamiltonian of the chain has the form: 



-4&W(t) 



E P< " (a (1) (t)e" 4 " (1)(t) + E (WPi + ftlft) + 

— 1 o V / — 1 1 



i=l,2 



i=l,2 



+ 



+ e" 



2(a«(t))V 4b(1) « + 72« 2 e 46(1) W + 4a«(t)6«(t) - d«(*)J E ^ + 36 «7 + 

V3, 



-4&«(t) 



7 ~°"37 

y\ 



7 2 - |7^"3 



[hi + f 1/2) 2 



+ 



7 



h a 3 



i=l,2 

2 ^r 



§2/! + 4 1/2) 2 



(35) 



The t-dependence of the kinetic term can be interpreted as a t-dependent mass 21 



Linear terms in momenta are known to describe for example the coupling of charged 
particles with gauge potentials and therefore have not to be discarded a priori. However 
the terms linearly dependent on momenta drop out for a particular relation: 



a(t) = -6(t)exp(46(t)). 



(36) 



We remark that, in the case of the factorization of h^ M \y,t) = r](t)h^ M \y), TDSE 
reduces effectively to a quasi- stationary problem, because by a suitable reparametrization 
of time t — > t = J r](t)dt the problem becomes stationary. The corresponding constraint 
leads for M = 0, 1, 2 again to (^) and to a nonlinear differential equation for b(t) : 



V(t) +66 2 (t) + 72a 2 (e' 



-86(t) 



-1} 



0. 



(37) 



The general solution of this equation involves elliptic integrals in the relation between t 
and b. The function T](t) becomes rj(t) = exp(— 46(i)). 

The construction of invariant operators R, which satisfy the equation: 



dR 

— +i[H($,t),R]=0 



(38) 



is an important aspect of the investigation of time-dependent systems [25]. In our frame- 
work from the intertwining relation (|30| ) the invariant operator exists and can be expressed 
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in terms of ft< M ' and U^: 

Ri M) (t) s uM{y,t)ht M \m {M) r l m= m 

where the last term is usually referred as gauge term. 

The invariant operator is hermitian because the intertwining operator U(y, t) is uni- 
tary. From the Eqs.(p£f), ( P5| ) and (^) it is straightforward to obtain the explicit expres- 



sion for the chain of invariants of this model. In particular, one can notice that have 
still the structure similar to the Calogero Hamiltonians ([34]), (|35D though some terms are 
missing. 

In general, one can argue from the similarity fl39D that the spectrum of R( M ^ is the 
same as the spectrum of and therefore time-independent ^7|. The operators R^ M ^ 
provide an additional exactly solvable (matrix and scalar) models with explicit time- 
dependent potentials but with time-independent spectra. Their eigenfunctions depend 
parametrically on time via U^ M \y, t) applied to the stationary eigenfunctions of hS M \ Let 
us remark that invariant operators R(t) admit a quasi-factorization like (|23j) in Section 
3 with suitable (transformed by U(y,t)) components of supercharge, but h(y,t) do not 
because of the gauge term. 



5. Conclusions 



Given for granted the usefulness of exactly (and quasi-exact ly) solvable models we 
would like to point out that our contribution has been to construct explicitly few models of 
such a kind with a discrete spectrum: among them the exactly solvable 3-particle (matrix 
and scalar) non-stationary Calogero models and quasi-exactly-solvable iV-particle matrix 
stationary models. An extension of the method of Section 4 to the iV-body Calogero 
model described in Section 3 leads to time-dependent quasi-exactly-solvable matrix mod- 
els. Since it is not usual to find exactly solvable or quasi-exactly solvable time-dependent 
problems, specially in a context of the many-body systems, our results support the pro- 
gram to investigate further time-dependent generalizations of stationary solvable models, 
like mentioned in the Introduction 0, ||, [|J, and quasi-exactly solvable matrix mod- 
els |2S| . In particular, one can focus attention on the dynamical algebras of these models 
p9| to construct Ermakov-Lewis invariant operators (|26| and references therein). A less 
straightforward task will be to modify the model in such a way as to allow for coexistence 
IpOfl of a continuum and a discrete spectrum describing scattering and bound states. 
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